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ABSTRACT 

Observations show the ubiquitous presence of propagating magnetohydrodynamic (MHD) kink waves in 
the solar atmosphere. Waves and flows are often observed simultaneously. Due to plasma inhomogeneity in 
the perpendicular direction to the magnetic field, kink waves are spatially damped by resonant absorption. 
The presence of flow may affect the wave spatial damping. Here, we investigate the effect of longitudinal 
background flow on the propagation and spatial damping of resonant kink waves in transversely nonuniform 
magnetic flux tubes. We combine approximate analytical theory with numerical investigation. The analytical 
theory uses the thin tube (TT) and thin boundary (TB) approximations to obtain expressions for the wavelength 
and the damping length. Numerically, we verify the previously obtained analytical expressions by means of 
the full solution of the resistive MHD eigenvalue problem beyond the TT and TB approximations. We find 
that the backward and forward propagating waves have different wavelengths and are damped on length scales 
that are inversely proportional to the frequency as in the static case. However, the factor of proportionality 
depends on the characteristics of the flow, so that the damping length differs from its static analogue. For slow, 
sub-Alfvenic flows the backward propagating wave gets damped on a shorter length scale than in the absence 
of flow, while for the forward propagating wave the damping length is longer. The different properties of the 
waves depending on their direction of propagation with respect to the background flow may be detected by the 
observations and may be relevant for seismological applications. 

Subject headings: Sun: oscillations — Sun: corona — Sun: atmosphere — Magnetohydrodynamics (MHD) 
— Waves 



1. INTRODUCTION 

The presence of ubiquitous small-amplitude propagating 
magnetohydrodynamic (MHD) waves in magnetic waveg- 
uides of the solar corona was first obser ved with the Coronal 
Multi-channel Polarimeter (C oMP) (see lTomczvk et alj2007t 
Tomczvk & Mclntoshll2009t) . The observations have inspired 
a number of recent theoretical works in which the proper- 
ties of propa gating waves are studied. Based on MHD wave 
theory (e .g.. lErdelvi & Fedunl 120071 IVan Do orsselae re et al.l 
I2008alfb1 : IVerth et al.l 12010b the observations have been in- 
terpreted as propagating kink waves, i.e., transverse MHD 
waves with mixed fast and Alfvenic properties whose 
dominating restoring force is magnetic tensi on (see, e.g., 
lEdwin & Roberts! 119831: IGoossens et al.l 12009ft . In particu- 
lar, |^n^oorsseiaereietaU(|2008b) showed that field-aligned 
density enhancements in the corona act as natural waveg- 
uides for MHD waves. Apart from the CoMP observations 
in coronal loops, propagating kink waves have also been 
observed in chromosp heric spicules (e.g., iDe Pontieu et al.l 
120071: iHe et aljj 2009a b) and i n thin threads of solar promi- 
nences (e.g.. iLin et al.l l2007, 2009|). In the same man- 
ner as standing kink MHD waves are damped i n time 
by resonant absorption (s ee, e.g., IGoossens et al.l 120021: 
iRuderman & Ro berts 2002), propagating kink MHD waves 
are damped in space due to naturally occurring plasma in- 
homo geneity in the directi o n transverse to the th e magnetic 
field dTerradas et aLl2 010b). Terrad as et al.l (1201 Obi) , hereafter 
TGV, obtained the important result that the damping length by 
resonant absorption is inversely proportional to the wave fre- 
quency. This means that high-frequency waves are damped in 

roberto.soler@wis.kuleuven.be 



shorte r length scales than low-frequency waves. IVerth et ail 
(2010) showed that the analytical theory of propagating res- 
onant kink waves developed by TGV is fully consistent with 
the CoMP observations . For the time-dependent, driven prob- 
lem lPascoe et al] (120101) studied numerically the spatial damp- 
ing of resonant kink waves and obtained equivale nt results to 
those a nalytically predicted by TGV. Importantly, Sol er et aT] 
d2011al) have shown that the results of TGV still hold when 
the plasma is partially ionized, so that the theory can be ap- 
plied to kink waves propagatin g in the chromosph ere and in 
prominences as well. Recently. ISoler et aD (1201 lbl) extended 
the results of TGV by taking into account for the first time 
density variation both transversely and along the magnetic 
field direction. 

Field-aligned flows are also ubiquitous in magnetic 
structures in the solar atmosphere (see t he observa- 
tional reports by, e.g . . iBrekke et al.l 119971 IZirker et al.l 
19981: IWinebarger et all 120011 120021 lOkamoto et al.l 120071: 
Chae et al.l 120081: lOfman & Wang|l2008l) . Typically, the ob- 



served flow velocities are smaller than 10% of the plasma 
Alfven speed. Faster flows of the order of the Alfven speed 
are much less frequent and are related to energeti c events as, 
e.g., flares and coronal mass ejections (see, e.g, llnnes et aTl 
2003). The investigation of the effect of flow on the proper- 
ties of the waves is therefore of evident interest. For example, 
the influence of flow and implications for MHD seismology of 
standing ki nk waves i n coro nal loops have been recen tly dis- 
cussed by Rude rmanl (120101) andlTerradas et al.| (l201 ll) . In the 
case of prominences, ISoler & Goossensl (1201 ll) investigated 
standing kink waves in coronal flux tubes partially filled with 
flowing threads of prominence material. Nevertheless, none 
of these works took damping into account. Also for stand- 
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ing waves, iTerradas et al.l (l2010al) studied temporal damping 
by resonant absorption in the presence of flow and found cor- 
rections to the dam ping time due to the flow with respect to 
the static case dGoossens et al.ll2002h . However, in the case 
of spatial damping of propagating kink waves the effect of 
flow has not been investigated. TGV did not include flow in 
their study. To our knowledge there is no work in the lit- 
erature that studies in detail the influence of flow on reso- 
nantly damped propagating waves in solar magnetic waveg- 
uides. Existing investigat i ons of surface waves in the solar 
wind as, e.g.. lEvans et al.l (120091) did not perform a rigorous 
treatment of the process of resonant absorption and consid- 
ered very simplified expressions for the damping length of the 
waves. Thus, a detailed investigation of propagating resonant 
MHD waves in the presence of flows in needed. 

Here we investigate the effect of flow on the spatial damp- 
ing of resonant kink waves in transversely nonuniform solar 
waveguides. We attack the problem both analytically and nu- 
merically. The analytical theory uses the thin tube (TT) and 
thin boundary (TB) approximations to obtain expressions for 
the wavelength and the damping length. We determine the 
influence of flow and compare our expressions with those ob- 
tained by TGV in the static case. Later, we use numerical 
methods to study the propagation and spatial damping of kink 
waves beyond the TT and TB approximations. The full nu- 
merical computations enable us to test the validity of the an- 
alytical expressions. Finally, we discuss the implications of 
our results for MHD seismology. 

2. MODEL AND GOVERNING EQUATIONS 

The equilibrium configuration is a straight cylindrical mag- 
netic flux tube of radius R embedded in a magnetized plasma 
environment. For convenience, we use cylindrical coordi- 
nates, namely r, <p, and z for the radial, azimuthal, and longitu- 
dinal coordinates, respectively. The z-direction is set along the 
axis of the cylinder. We use the B = approximation, where B 
is the ratio of gas pressure to magnetic pressure. The B — ap- 
proximation enables us to arbitrarily choose the density pro- 
file. In what follows, subscripts i and e refer to the internal 
and external plasmas, respectively. For example, we denote 
by pi and p e the internal and external densities, respectively. 
Both of these quantities are constants. There is a nonuniform 
transitional layer in the transverse direction that continuously 
connects the internal density to the external density. The layer 
has a thickness / and covers the interval R-l/2 < r < R + 1/2. 
The equilibrium magnetic field is straight, B = Be,, with B 
constant. We assume a background flow along the magnetic 
field direction, U = Ue z . We denote by E/j and U e the internal 
and external flow velocities, respectively, which are constants. 
We take U\ and U e as positive quantities. As for the density, 
we allow the flow velocity to change continuously in the ra- 
dial direction from its internal to its external values within a 
transitional layer. The transitional layer for the flow velocity 
extends in the interval R - I* /2 < r < R + I* /2, with I* the 
thickness of the transition. 

Linear, ideal MHD waves propagating in our model are 
governed by the following set of equations, 

p|tt +U- Vv + v- VU|=i(Vxb)xB, (1) 
\dt )n 

^ - Vx(Uxb) = Vx(vxB), (2) 
ot 

where p is the plasma density, v is the velocity perturbation, 



b is the magnetic field perturbation, and p is the magnetic 
permittivity. As the equilibrium is uniform in both cp- and z- 
directions and we consider waves propagating along the tube 
with a fixed frequency, we write all perturbations proportional 
to exp (im<p + ik z z — iu>t), where m is the azimuthal wavenum- 
ber, k z is the longitudinal wavenumber, and u> is the wave 
angular frequency. Due to the presence of a transverse in- 
homogeneous transitional layer, wave modes with m + are 
spatially damped due to resonant absorption. Here we are in- 
terested in kink waves, which are described by m = I. Kink 
waves have mixed Alfvenic and fast MHD properties. They 
are the only wave modes that can displace the magnetic cylin- 
der axis and s o produce transverse motions of the whole flux 
tube (see, e.g jEdwin & Roberts! 1983UGoossens et al.ll2009l) . 
As a result of the process of resonant absorption, transverse 
kink motions of the flux tube are damped and azimuthal mo- 
tions within the transitional layer are amplified as the wave 
propagates along the magnetic cylinder. 

For real w, resonant damping causes k z to be complex, 
k z = k z n + ik z \, with k z & and k z \ the real and imaginary 
parts of k z , respectively. For fixed and positive u>, the di- 
rection of wave propagation is determined by the sign of 
k Z R. For £.r > the wave propagates towards the posi- 
tive z-direction (forward waves), whereas for k z & < the 
wave propagates towards the negative z-direction (backward 
waves). In the absence of flow both directions of propa- 
gation are equivalent. In the presence of flow forward and 
backward waves have not the same properties and both di- 
recti ons of propagation must be taken into account (see, 
e.g., [N akariakov & Roberts] 1 19951: iTerra-Homem et al.l 120031 
ISoler et al.ll2008t IVasheghani Farahani et al.ll2009l) . Regard- 
ing the imaginary part of k z , resonant absorption spatially 
damps the wave, so we expect k z i > 0. However, strong flows 
may trigger the Kelvin-Helmholtz I nstability (KHI) (see, e.g., 
Chandrase kharll96TtlDrazin & Reidll 19811) . causing modes to 
be amplified in z, i.e., k z \ < 0. From the real and imaginary 
parts of k z we compute the wavelength, A, and the damping 
length, Ld, as 

2n 1 
-4=7-, L» = —. (3) 

k;R k zl 

3. ANALYTICAL INVESTIGATION 

To study analytically the effect of flow on the resonantly 
damped kink waves we use the TT and TB approximations. 
In the TT approximation we restrict ourselves to waves with 
A/R » 1. In terms of frequency, the TT approximation is 
equivalent to the low-frequency approximation, i.e., wta <k 
1, with ta = R/va the Alfven travel time and va = Bj -\fpp 
the Alfven velocity. The TB approximation is used here to in- 
clude the effect of resonant absorption in the inhomogeneous 
layer, and is valid for l/R <K 1 and I* /R «: 1. In the TB ap- 
proximation, the jump of the perturbations across the inhomo- 
geneous layer is assumed to be the same as their jump across 
the resonant layer. Th e expressions f or the j ump conditions 
can be found in, e.g.. ISakurai et all (119911) : iGoossens et al.l 
d 1 9951): iTir rv & Goossens (1996) for the sta tic case, and in 
IGoossens et alj dl992F lErdelvi et aj] (119951) for the station- 
ary case. Then, the connection formulae at the Alfven res- 
onance are used as jump conditions for the perturbations at 
the tube boundary (see extensive deta i ls about the method i n 
IGoossens et~aT]|2006l: lGoossensll20()8t IGoossens et~alll2011h . 
The dispersion relation for kink and fluting MHD waves in 
the TT and TB approximations is (see, e.g., IGoossens et al.l 
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p(r A )|A A 



Pi(of-Op e (o2- W y ( 4) 



where Q = a> - Uk z is the Doppler-shifted frequency, a> A = 
k\v 2 K is the square of the Alfven frequency, r A is the Alfven 
resonance position, and 



— [Q 2 



(5) 



The term on the right-hand side of Equation © contains the 
effect of resonant absorption. Note that the dependence on 
m of the dispersion relation is only present in this term. This 
means that in the TT approximation the effect of the azimuthal 
wavenumber, m, is only felt in the damping of the wave, not 
in its propagation. 

In the case of temporal damping of standing waves, i.e., 
real k z and complex co, the sol utions to Equat ion <j4j have 
been investigated in detail by Goossens et al. (1992) and 
iTerradas et al.l(l2010al) . In the case of spatial damping of prop- 
agating waves, i.e., real u> and complex k z , Equation has 
been explored by TGV in the absence of flow. Here our pur- 
pose is to investigate the effect of flow on the results obtained 
by TGV in the static case. 

3.1. No resonant damping 

In the absence of resonant damping, i.e., for l/R = I* /R = 
0, Equation © becomes 

ofa) (fa) =0. (6) 

Equation © is independent of m. In the static case, C/j = 
U e — and O; = Q. e — u>. The solution to Equation © is 



k z = ±— = ±ko, 

Vk 



with 



v k = 



/ 2 2 \l/2 



(7) 



(8) 



Pi +Pe 

the kink velocity. In Equation (|7]l the + sign stands for the 
forward wave and the - sign for the backward wave. In the 
absence of flow both forward and backward waves are equiv- 
alent and have the same wavelength, A = 2njk^. 

In the presence of flow, we rewrite Equation (|6) as a second- 
order polynomial in k z , namely 



( v k ~ v ke) + 2v cm w k, - u? = 0, 
where v cm and vke are defined as 

( Pl Uf + p e U 2 e ^ 2 



p-JJi +p e U e 



vke = 



Pi +Pe 



(9) 



(10) 



Vom is the center-of-mass velocity, vke is the velocity asso- 
ciated to the kinetic energy of the flow. Note that, similarly, 
the kink velocity, Vk, is the velocity associated to the energy 
of the magnetic field. The two solutions to Equation (O for 
v K E * v k are 

1/2 



k, = - 



-to+ko 



1 - 



PiPe (Ui - u e y 



(Pi + Pe) 
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where the + and - signs in front of the second term stand for 
forward and backward waves, respectively. If C/j — t/ e = 0, 
Equation (fTTT > simply reduces to Equation (O. We clearly see 
in Equation (ITTb that the equivalence between both directions 
of propagation is broken by the flow. For vke < v k forward 
and backward waves propagate in opposite directions. For 
vke > ^ both waves propagate in the same direction, i.e., they 
both are forward waves in practice because the flow is strong 
enough to force the backward wave to reverse its direction of 
propagation. In the particular case vke = Vk, the solution to 
Equation (O is 

K = (12) 

2v cm 

which corresponds to the forward wave, while the backward 
wave does not propagate in the static reference frame. For 
U e — the condition Vke - v k is equivalent to 



Uf = 2v 



Ai' 



(13) 

When the argument of the square root in Equation dm is 
negati ve, k 7 becomes complex. This is the classic al KHI (see, 
e.g.. iChandrasekhari [19611 iDrazin & Reidlll981l) . Then, the 
two solutions correspond to a spatially damped mode and a 
spatially amplified mode, respectively. The KHI appears for a 
critical velocity shear, AU = C/j - U e , defined as 



(At/) 2 > 



(p +p e y 

PiPe 
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v,, = v v 



Again, in the reference frame where U e 
be rewritten as 

Pi 



Uf > 2 1 + 



Pe 



(14) 

Equation ( TT4l > can 
(15) 



From Equations ( fT~3T > and ( fT~5T > we see that the KHI requires a 
faster flow velocity than the one needed to reverse the propa- 
gation of the backward wave. 

Equivalently to the case without flow (see Equation (0), we 
can rewrite Equation (ITTb as 



v kf 



(16) 



where we have introduced the effective kink velocity modified 
by the flow, vi, defined by 



v k 



k kf 



PiPe (Ui - u e y 



(Pi +Pe) 



1/2) 



(17) 



Note that the effective kink velocity is different for forward (+ 
sign) and backward (- sign) waves. This means that the phase 
speed of the waves depends on their direction of propagation. 
Also note that becomes complex beyond the critical veloc- 
ity shear for the KHI. For slow, sub-Alfvenic flows we may 
drop the quadratic terms in Ui and U e from Equation (TTTb to 
obtain a first-order approximation for vi, namely 

v£f « ±v k + v cm . (18) 

In the absence of flow, Equation ( TTTb reduces to vj^ = ±Vk. 
When the effective kink velocity is equal to the external 
Alfven velocity, the wave becomes leaky in the external 
medium. In terms of k z , waves are leaky for wavenumbers 
larger than 

k z = ±— = ±k h . (19) 

v Ae 
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The forward wave becomes leaky for much slower flow ve- 
locities than the backward wave. By using Equation ( TT~8T > and 
taking U e = 0, we find that the forward wave becomes leaky 
for 
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Pi J 
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(20) 



Again, we use the approximation of Equation ( TT~8T > in Equa- 
tion ([Tol l to obtain a first-order approximation of the wave- 
length as 



A*/lo ±1 + 



(21) 



with /lo = 2n/ko. We consider the particular case t/ e = 
and use the dimensionless notation of TGV to rewrite Equa- 
tion (fJTJ as 



A 
R 



2n 



2( 1 



\f +1/ 



±1 + 



2(f +1) 



f/i 



with 



f. _ Pi. /= — 

Pe ' V A i ' 



VAi 



(22) 



(23) 



the density contrast, the dimensionless frequency, and the di- 
mensionless flow velocity, respectively. 

We plot in Figure 03a) \A\/R versus U\ computed from the 
full Equation (fTlT i for the particular case U e = 0. As predicted 
by Equation d!31 l. the backward wave reverts its direction of 
propagation for U- t = V2. For the set of parameters used in 
Figure [TJ the forward wave becomes leaky for a flow velocity 
slightly sub-Alfvenic (Equation (T20li). When the threshold ve- 
locity of the KHI is reached (Equation ifTBT)). both forward and 
backward waves merge. We compare the full result with the 
approximation for slow flows given by Equation (1221 (see the 
symbols in Fig. Q2 a )X an d obtain a good agreement for sub- 
Alfvenic flows, i.e., Ui < 1. On the other hand, Figure [Tib) 
displays k z \R versus Ui. In the absence of resonant damping, 
the imaginary part of k z is zero for flow velocities slower than 
the critical velocity shear of the KHI. For larger velocities, 
one damped solution and one overstable solution are present. 

3.2. Effect of resonant damping 

Here we incorporate the damping due to resonant absorp- 
tion. We take into account the full expression of the disper- 
sion relation (Equation ©). We write k z = k zR + ikg in Equa- 
tion dU and consider weak damping, so we neglect terms of 
O In addition, we implicitly assume that the flow veloci- 
ties are slower than the critical velocity of the KHI. After long 
but straightforward analytical manipulations, we obtain from 
Equation (0]i the expression for the ratio k z i/k zR , namely 



1 ~ ^Ai) 2 



(24) 



2 r A pi + p e p (r A ) | A A | O)(eo- Lo cm ) 

with o) cm = kzRV cm . In the absence of flows, Equation d24i > 
reduces to Equation (8) of TGV. Due to the different values of 
Q; and a> cm for forward and backward waves, Equation (l24l 
predicts that both waves have different damping ratios. The 
expression for the ratio k z j/k Z R is more complicated in the 
presence of flows compared to the static case of TGV. Let us 
find a more simple expression for k z i/k zR in the limit of slow, 
sub-Alfvenic flows. First, we evaluate p (r A ) |A| A to explicitly 
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Figure 1. (a) \ A\ /R and (b) k z \R versus U\ corresponding to the forward (solid 
line) and backward (dashed line) kink waves in the absence of resonant damp- 
ing. The vertical dotted lines in both panels correspond to the different critical 
flow velocities indicated in the text. The horizontal dotted line in panel (a) is 
the wavelength for V\ = 0, with Aq = 2n/ko. The sym bols correspond to the 
approximation for slow flows given in Equation i22\ . In panel (b) the hori- 
zontal dotted line denotes k z \ = 0. In this plot, / = 0.1, £ = 3, and U c = 0. 



take into account the radial variation of the density and the 
flow velocity at the resonance position. From Equation (|5j 
we get 



p(r A ) |A A | 



Q 2 (r A ) 



- 2p (r A ) Q. (r A ) k zR 



'dU\ 



(25) 

where we have used the resonant condition Q. 2 (r A ) = 
k 2 R v 2 A (r A ). The quantities present in Equation (f25l have to 
be evaluated at the resonance position, r A . In the TB approx- 
imation it is reasonable to assume r A m R. There are two 
terms in the right-hand side of Equation (l25T l. The first term 
is due to the variation of density and the second term is due 
to the variation of flow velocity. The direction of wave prop- 
agation is also important in Equation ( |25T l because of the sign 
of k Z R. In the absence of flow, Equation d25l > simplifies to 

p(r A )|A A |=^(g) rA . 

For smooth profiles, the derivatives of the density and the 
flow velocity profiles at the resonance position can be cast as 



dp\ „, T n Pi ~Pe 
dr) 4 I ' 



«L\ ,T«-^ (26) 
dr L 4 I* V ' 



with T a factor that depends on the form of the trans- 
verse profile. For ex ample, T = 4/jt 2 for a linear profile 
dGoossens et al.l 120021) and T = 2 In for a sinusoidal profile 
dRuderman & Roberts 2002). For simplicity we assume the 
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same profile for both density and flow velocity, but we keep 
I + /*. As values of density and flow velocity at the resonance 
position we take 



pO*a) 



Pi +Pe 



U(r A ) 



Ui + U e 



(27) 



We use in Equation d25l the expressions given in Equa- 
tions d26l and d2Tb to obtain 



p(r A )\A A \=T 
x 1 



X Pi~Pe 2 

4 Z 



(r A ) 



fc;R Pi + Pe 
^ ('"a) Pi - Pe 



(t/i 



(28) 



For our following analysis, we assume that the second term 
within the absolute value in Equation ( 1281 is smaller than 
one, so that we can drop the absolute value sign. For slow 
flows, this is equivalent to assume I* > I. In the limit of 
slow, sub-Alfvenic flows, we neglect the quadratic terms in 
the flow velocities. In addition, we write k zR * co/v^ and use 
the first-order approximation for v^ f given in Equation (list . 
Equation ( 1281 can then be rewritten as 



t MAI ar n A ~ Pe 2 

p (r A ) | A A | « T- : — co 



X 1 + 



4 / 
Ui + U e 



1 + 



Pi +p e Ui-U e I 



Pi - Pe t/i + U e I* 



■(29) 



Next, we use Equation d29l in Equation ( l24l and perform 
a first-order expansion in the flow velocities. Equation (l24l 
becomes 



1 m I pi -p e j i', 



In T R p\ + p & 



1-4 



Pi +Pe PiUi~p e U e 



Pi -PePiUi +p e U e 



Ui-U e 



1 + 



Pi + p e t/j - I/e / 
Pi -Pe t/i + t/e /* 



(30) 



As in previous expressions, the + and - signs in Equation (f30b 
stand for forward and backward waves, respectively. In the 
particular case U e - 0, Equation (1301 simplifies to 



1 m I p; - p e 



27T y /? p; + p e 

v k \ Pi + Pe Pi - Pe Pi - Pe <* 



4p; Pi + p e / 



(3D 



We use again the approximation Ic-r * co/vtf and compute 
from Equation (1311 the damping length, Ld = 1 /&zi, as 



T Rp\+p e v k 
2n— 

m I pi - p e U) 



j ± tA / 3pi + Pe 



Vk \ p 



Pi + Pe j 
Pi - Pe /* 



(32) 

Finally, we express Equation ( [321 using the dimensionless 
quantities defined in Equation (|23l , namely 



with 



2£ 1 



\f +1/ 



1 ± t/;. 



f + 1 /3£+ 1 f +1 / 



2f \ f - 1 f - 1 /" 



^ E m / £ - 1' 



(33) 



(34) 



Equation d33l is the key equation of this investigation and 
contains basic properties on the spatial damping of propagat- 
ing kink MHD waves. Several important results can be ex- 
tracted from Equation (l33l . Equation (l33l l predicts that back- 
ward and forward propagating waves are damped on length 
scales that are inversely proportional to the frequency, /. This 
is the same dependence found by TGV in the static case. 
However, the factor of proportionality depends on the char- 
acteristics of the flow and the density contrast, so that the 
damping length differs from its static analogue. As for the 
wavelength (see Equation (l22l >). the damping length for for- 
ward and backward waves is different. 

To shed more light on this result, let us consider the case 
I* » I. Equation d33l reduces to 



27Tft 



2<T 1 



1 ± Ui- 



£ + j3f+j 
' 2f £-1 



(35) 



According to Equation (1351 1 the backward propagating wave 
(- sign) gets damped on a shorter length scale than in the 
absence of flow, while for the forward propagating wave (+ 
sign) the damping length is longer. For I * /* the damping 
length for the forward wave remains longer than that of the 
backward wave, but we must note that for /* <s / the situation 
may be the opposite. From Equation (l33l l we can assess the 
relation between / and /* for which the two terms multiplying 
the flow velocity cancel each other, namely 



3^+1 



(36) 



Thus, the backward wave damping length becomes longer 
than that of the forward wave for I* smaller than the value 
given in Equation d36l . Strictly, Equation d33l is not valid in 
the limit I* <K /, because we have to properly take into ac- 
count the absolute value in Equation ( 1281 . For slow flows and 
I* <k /, Equation d28l can be approximated as 



( \ I A I ar n Pi +Pe 2 Uj-Ue 

p (r A ) |A A | ~ T- — - — co 

4 I* Vk 



(37) 



Importantly, we find that Equation (T3~7l is independent of I and 
the same expression holds for forward and backward waves. 
Now, from Equation (T3~7l it is straightforward to obtain that 
Ld ~ 1 //* for both forwa rd and backward waves in the limit 
I* <K /. As discussed by Terrada s et al.l d2010a) in the case 
of temporal damping of standing waves, multiple resonances 
may occur within the inhomog eneous transitional layer in the 
limit I* <K I (see Fig. 5 of Terradas et al. 2010a). In such a 
case, the total damping rate is the sum of the contributions 
from each resonance. We do not study this peculiar situation 
which takes place for very small, probably not rea listic I*. In- 
stead, we refer the reader to Terrada s~et al.l(B0 10a) for details. 

4. NUMERICAL COMPUTATIONS 

Here, we verify the validity of the analytical expressions 
derived in Section [3] in the TT and TB approximations. To 
do so, we numerically solve the full eigen value problem by 
means of the PDE2D code (ISewellll200l . The numerical 
scheme is similar to that used by TGV. The code implements 
a method based on finite elements to numerically integrate 
Equations ([1} and (f2]i in the radial direction from the cylinder 
axis, r = 0, to the edge of the numerical domain, r = r max . 
The boundary conditions at r — are set according the sym- 
metry arguments, while we impose all perturbations to vanish 
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at r — r max . In order to obtain a good convergence of the 
solutions, r max is located far enough from the magnetic tube 
to avoid numerical errors. We take r max = 100/?. We use a 
nonuniform grid with a large density of grid points within the 
inhomogeneous layer in order to correctly describe the small 
spatial scales of the eigenfunctions due to the Alfven reso- 
nance. To avoid the singularity of the ideal MHD equations 
at the resonance position, we add to the induction equation 
(Equation (O) a resistive term, i.e., rjV 2 b, with 77 the mag- 
netic diffusivity. The output of the code is the complex eigen- 
functions and their corresponding eigenvalues. In the limit 
of large R m , with R m = v^R/t] the magnetic Reynolds num- 
ber, the eigenvalues are independent of R m . In such a case, 
wave damping is due to resonant absorption exclusively. In 
our computations, we have considered a sufficiently large R m 
and have checked that the eigenvalues are indeed independent 
of R m . We typically take R m * 10 7 in the computations. 

The PDE2D code solves the eigenvalue problem for the 
temporal damping, i.e., for complex u provided a fixed and 
real k z . As in TGV, we need to convert the results from tem- 
poral damping to spatial damping. We use Equation (40) of 
TGV to perform the conversion from temporal damping (com- 
plex a) and real k z ) to spatial damping (real a> and complex k z ). 
From Equation (40) of TGV it is straightforward to obtain the 
relation between the imaginary parts of a> and k z , namely 

where a> R and <x>\ are the real and imaginary parts of the fre- 
quency in the temporal damping case. In the spatial damp- 
ing case cj = wr. The factor within the parenthesis is the 
group velocity. Note that the flow does not explicitly appear 
in Equation d38l l. but its effect is contained in the value of 
group velocity numerically computed. 

First of all, we test the numerical code by considering the 
case without flow. In this case we fully recover the results 
of TGV. Thus, we are confident that the code works properly. 
Hereon we incorporate the effect of flow. Regarding the wave- 
length, our numerical results indicate that the presence of the 
transitional layer has a small impact on the value of the wave- 
length. The wavelength in the case without resonant damping 
(see Fig. [Tf a)) is a very good approximation to the wavelength 
in the resonant case. The approximation for slow flows given 
in Equation (1221 also holds in the case with damping. There- 
fore, our following analysis is focused on the behavior of Lq 
in the presence of flows. The numerical results are used to test 
the approximations behind Equation (l33l l. i.e., the TT and TB 
approximations, and the assumption of slow flows. 

In Figure [2] we display Lu/R as a function of the dimen- 
sionless frequency, /, for U\ - 0.1 (the rest of parameters 
are indicated in the caption of the Fig. |2). As in the static 
case of TGV, the higher the frequency, the shorter the damp- 
ing length by resonant absorption. We obtain the analytically 
predicted result that forward and backward waves have differ- 
ent damping lengths. In this example, the damping length of 
the forward wave is longer than that of the backward wave. 
The equivalent damping length in the absence of flow is in 
between both values (the dotted line in Fig. |2}. We compare 
the numerical results with those in the TT approximation and 
for slow flows (Equation (l33l). The TT approximation applies 
when / <k 1. In Figure |2j a) we plot the results for / < 0.1. 
An excellent agreement is found between numerical and an- 
alytical solutions for both forward and backward waves. The 




0.0 0.2 0.4 0.6 0.8 1.0 

J 



Figure 2. (a) Ratio of the damping length to the radius, L^/R, versus the 
dimensionless frequency, /, corresponding to the forward (solid line) and 
backward (dashed line) kink waves for U\ = 0.1, U e = 0, f = 3, and l/R = 
I* /R = 0.1. The sym bols correspond to the approximation for slow flows 
given in Equation )33t . The dotted line is the result in the absence of flow 
(see TGV). (b) Same as panel (a) but for larger /. 

agreement between approximate and numerical results is also 
remarkably good even when the condition / « 1 of the TT 
approximation is not strictly fulfilled. This can be seen in Fig- 
ure|2jb), where results are plotted for / < 1. 

Next we determine the influence of / and I* , and test the 
TB approximation used to derive the analytical expressions. 
First we consider the case 1 = 1*. Figure [3} a) shows the 
dependence of the damping length on l/R. As for the fre- 
quency, Lu/R decreases as l/R increases. The dependence 
of the damping length on l/R for both forward and backward 
waves is the same. As before, a very good agreement be- 
tween Equation ( T33T > and the numerical result is found even 
when the thickness of the transitional layer departs from the 
limit l/R <k 1. This means that the TB approximation is suf- 
ficiently accurate when the condition l/R <k 1 is slightly re- 
laxed. This result enables us to confidently use Equation < f33T > 
beyond the condition of validity of the TB approximation. Al- 
ternately, in Figure|3lb) we keep l/R constant and vary I* /R, 
so that the spatial scales for the variation of density and flow 
velocity are different. We recover the analytically predicted 
result that for I* » / the results are independent of I* /R. For 
I* <sc I Equation ( f33l > does not correctly describe the damping 
length of the forward wave. As discussed at the end of Sec- 
tion [372J Ld ~ 1//* in the limit I* <k /, and so Lu/R increases 
as l*/R 0. 

Finally, we assess the effect of the flow velocity. This is 
done in Figure [4] Again, the linear approximation (Equa- 
tion (l33l ) is quite accurate and agrees well with the full 
numerical results for t/i < 0.1. As expected, the differ- 



RESONANT KINK WAVES WITH FLOW 



7 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 

l/R U- 



( b ! 1500 




Figure 3. (a) Ratio of the damping length to the radius, L&/R, versus l/R 
for the forward (solid line) and backward (dashed line) kink waves in the 
case I* /R = l/R. The dotted line is the result in the absence of flow, (b) 
Dependence on I* /R for l/R = 0.15. The vertical dotted line denotes the 
value of I* /R fo r wh ich both forward and backward waves have the same 
Ld/R (Eq uati on )36» . In both panels the symbols are the approximation of 
Equation {33). In these plots, U; = 0.05, U e = 0,f = 0.1, and <f = 3. 



ence between the numerical results and the linear approxi- 
mation increases as the flow velocity gets faster. However, 
for Ui = 0.2 the relative difference between the full solu- 
tion and the approximation is only around 10% for the for- 
ward wave and 30% for the backward wave. This means 
that for rea listic flow velocities observed in coronal magneti c 
loops (e.g., Bre kke et alJl997tlWinebarger et al.ll20 01. 2002), 
Equation (l33l correctly describes the behavior of the damping 
length. 

In summary, in this Section we have confirmed that the ana- 
lytical expressions obtained in the TT and TB approximations 
and for slow, sub-Alfvenic flows are very accurate even when 
these expressions are used outside their domain of strict va- 
lidity. This result enables us to use Equation ( 1331 . i.e., the 
key equation of this investigation, when realistic values of 
frequency, flow velocity, and the rest of relevant parameters 
obtained from the observations are used. 

5. DISCUSSION 

Naturally, kink waves propagating in nonuniform magnetic 
flux tubes are spatially damped by resonant absorption. In the 
static case, TGV showed that the damping length is inversely 
proportional to the frequency. Here we have investigated ana- 
lytically and numerically the spatial damping of resonant kink 
waves in a transversely nonuniform magnetic waveguide in 
the presence of longitudinal background flow. Longitudinal 
flow breaks the equivalence between forward and backward 
propagating waves with respect to the flow direction. The 



Figure 4. Ratio of the damping length to the radius, L&/R, versus the flow 
velocity normalized to the internal Alfven velocity, U\, for the forward (solid 
line) and backward (dashed line) kink waves. The symbols are the linear 
approximation given in Equation l |33t . We have used l/R = I* /R = 0.1, 
r/ e = 0,/ = 0.1,andf = 3. 



wavelength and the damping length due to resonant absorp- 
tion are both affected by the flow. For sub-Alfvenic flows, the 
backward wavelength is shorter than that of the forward wave, 
and backward waves are damped in shorter length scales than 
forward waves. However, as in TGV we have found that the 
damping length of both forward and backward propagating 
waves is inversely proportional to the frequency. 

MHD seismology based on propagating waves has at- 
tracted limited attention and definitely less than its counter- 
part based on standing kink waves. Standing kink MHD 
waves are rare phenomena as they need a violent and en- 
ergeti c event such as a solar flare for their excitatio n (see, 
e.g.. lAschwanden et al.|[l999tlNakariakov et al.lll999l) . In the 
absence o f flow and for coronal loop standing oscillations 
(see, e.g.,[N akariakov & Ofman 2001; Goos sens et al.ll2002T: 
lArregui et al.ll2007l l2008UGoossens et al.ll2008l) . MHD seis- 
mology has been used to obtain information of the plasma 
physical conditions. Particularly, for a given set of parame- 
ters provided by the observations, i.e., period, damping time, 
and w avele ngth in the case of stan ding waves, Arreg ui et al.l 
(120071) and IGoossens et al.l d2008l) showed that the possible 
values of VAi, £, and l/R which are consistent with the the- 
ory form a one-dimensional curve in the three-dimensional 
parameter space. In principle, any point of th is curve can 
equally explain the obser vations. ISoler et all (120101) and 
Arregui & Ball esterl (1201 ll) showed that more constrained es- 
timations of v'Ai and l/R can be given in the case of promi- 
nence thread os cillations as the hmit C » 1 can be adopted. 
More recently, Arr egui & Asensio Ramosl (1201 ll) found that 
more accurate estimations of the par ameters are poss i ble by 
combining the analytical theory of IGoossens et al.l (120081) 
with statistical Baye sian analysis. In the presence of flows, 
iTerradas et al.l (|201 1) have recently explained also for stand- 
ing waves that the flow velocity can be estimated from the 
wave phase difference along the magnetic loop. 

On the contrary, propagating MHD wa ves are ubiq- 
uitous in the solar atmosphere (see, e.g., iTomczvk et al.l 
120071: ITomczvk & Mcinto sh 2009) and provide a huge reser- 
voir of possibilities for seismology. Some examples of 
MHD seismology bas e d on propagating waves are, e.g., 
IVan Do orsselae re et al.l (|2008b) using numerical simulations 
of guided MH D waves b y density enhancements in the so- 
lar corona, iLin et all (2009) using observations of kink waves 
in prominence threads, Vert h et aLl (120101) using resonantly 
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damped kink waves in coronal loops, and IVerth et all (1201 lb 
exploiting the properties of kink waves in chromospheric 
spicules. However, none of these works included flow in their 
analysis. Our theoretical results given in Equations (l22t and 
( 1331 ) have direct implications for MHD seismology based on 
propagating waves in a flowing medium, and could be used to 
infer information about the plasma properties. Therefore, the 
potential application of MHD seismology to the case of reso- 
nantly damped propagating kink waves in a flowing medium 
must be explored. 

In the presence of flow two waves with different wave- 
lengths and damping lengths but with the same frequency (or 
period) are simultaneously present. We denote as A + and Lu + 
the wavelength and damping length of the forward wave, re- 
spectively, and as X~ and Lp the equivalent quantities of the 
backward wave. Observationally, this means that it is possible 
to measure five quantities, namely A + , A~, Lq + , Lq , and the 
period, while the rest of parameters, i.e., VAi, U\, l/R, and ( are 
in principle unknown. In this analysis we assume l/R = I* /R 
for simplicity. We also take A~ as a positive quantity and so we 
perform the absolute value of Equation (l22l when the - sign 
is used. If observations can provide us with reliable values for 
the wavelengths and damping lengths of the two waves, then 
we can use our theoretical results (Equations (l22t and (1331 ) to 
obtain seismological estimations for the unknown quantities 
VAi, Uu l/R, and £ as 



v a = - 




P { 



2 

A + +A- 



t - 1 L D + + L D 
1 +2y 



\-2j 

with P = 2n/a> the period and 

A + — A Ld + + 



7 = 



A + + A- L D + - L D - 



(39) 

(40) 
(41) 
(42) 

(43) 



In the absence of flow, A + = A~ and Lq + = ■ Then, Equa- 
tio ns J39l and (PiTT i are equivalent to the expressions studied 
by Goosse ns et al.1 (120081) . and the density contrast (Equa- 
tion ( l42l i) becomes indeterminate. Thanks to the flow, the 
density contrast can be determined if reliable measures of the 
wavelengths and damping lengths of both forward and back- 
ward waves are available. 

The theoretical results of the present paper offer new and 
exciting opportunities for MHD seismology in plasma struc- 
tures with equilibrium flows. For the first time we have shown 
that an estimation of the density contrast, £ = pi/p e , is possi- 
ble. MHD seismology requires theory and observations. The 
required observations might not be available at present time. 
However, the seismological tool provided in Equations (|39l- 
( l42l i could be used in the future when the required observa- 
tions become available. 

The investigation performed in this paper may be improved 
in the future by incorporating additional physics in the MHD 
wave model. Effects that come to mind are the variation of 
the plasma paramet ers along the magnetic field direction as in 
ISoler et all (1201 lbl) and magnetic expansion and twist of the 



flux tube. These and other effects might be included in forth- 
coming investigations on propagating resonant kink waves. 
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